WAVES IN RADIAL GRAVITY USING MAGNETIC FLUID. 

D. R. Ohlsen 1 , J. E. Hart 2 , and P. D. Weidman 3 , 

’Program in Atmospheric and Oceanic Sciences, University of Colorado, Boulder, CO 80309-031 1, 

ohlsen@fluida.colorado.fidu 

2 Program in Atmospheric and Oceanic Sciences, University of Colorado, Boulder, CO 80309-03 11, 

hart@tack.colorado.edu 

department of Mechanical Engineering, University of Colorado, Boulder, CO 80309-0427, 

Patrick.Weidman@coIorado.edu 


INTRODUCTION 

Terrestrial laboratory experiments studying various 
fluid dynamical processes are constrained, by being in 
an Earth laboratory, to have a gravitational body force 
which is uniform and unidirectional. Therefore fluid 
free-surfaces are horizontal and flat. Such free surfaces 
must have a vertical solid boundary to keep the fluid 
from spreading horizontally along a gravitational po- 
tential surface. In the absence of terrestrial gravity, 
surface tension forms fluid masses into spherical balls 
without solid boundaries, as demonstrated on the 
Space Shuttle in the Drop Physics Module [1]. A 
fundamentally different problem is the behavior of flu- 
ids with a body force rather than a surface force that 
generates the spherical geometry. In many atmos- 
pheric, oceanic, or stellar fluid flows, the horizontal 
scale is large enough that sphericity is important in the 
dynamics. Further, fluids in spherical geometry can 
cover an entire domain without any sidewall effects, 
i.e. have truly periodic boundary conditions. We de- 
scribe novel spherical body-force laboratory experi- 
ments using ferrofluid. 

Ferrofluids are dilute suspensions of magnetic di- 
poles, for example magnetite particles of order 10 nm 
diameter, suspended in a carrier fluid [2]. A surfactant 
coating keeps the particles separate so that thermal 
Brownian motions in the fluid are sufficient to over- 
come both gravity and particle-particle attraction to 
keep the dipoles in suspension. For flows in which 
external magnetic field variations and the bulk fluid 
motions are slow compared to the time for the mag- 
netic fluid particles to rotate (~10 6 s), the fluid mag- 
netization, M , is parallel to the applied magnetic field, 
ff. If there is a gradient in H- \H\, then there will be 
a systematic body force owing to the slight but persis- 
tent correlation of field strength and pole sense. With 
this additional body force term, the momentum equa- 
tions become [2]: 

— + w V«+2Qxw = - — + g + -^^- + vV 2 u (1). 
dt p p 

The superparamagnetic ferrofluid response is com- 
plicated by the fact that M = M{H) unlike p which is 
independent of g for the gravitational case. As the ap- 
plied magnetic field, H , increases, more of the ferrofluid 
dipoles line up with the field resulting in a linear in- 
crease in ferrofluid magnetization, M , for small H but 


saturating at M s for large H. Figure 1 shows an exam- 
ple measurement of the field-averaged magnetization, 
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that shows this effect. 
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Figure 1. Field averaged magnetization vs. 

magnetic field for diluted water-based ferrofluid. 

The curve is a least-squares fit to a theoretical 

Langevin function prediction [2, 3], 

In the experiments described below, the magnetic 
fields are less than 400 Gauss so M is nearly propor- 
tional to H. We have made careful measurements of M 
vs. H for H < 500 Gauss and fit them empirically to 
polynomials in H to use in theoretical predictions of 
the body force, MVH. 

APPROACH 

The present study improves on the laboratory tech- 
nique introduced in [3]. Figure 2 shows a schematic of 
the experiment. The apparatus is cylindrical ly axi- 
symmetric. A cylindrical ceramic magnet is embedded 
in a smooth, solid, spherical PVC ball. The body 
forceon the water-based ferrofluiddepends only on the 
field strength and gradient from the interior magnet. 
The geopotential field (surfaces of constant | H |) and 
its gradient, the body force, were made nearly spherical 
by carefulchoice of magnet heigh t-to-diameter ratio and 
magnet size relative to the PVC ball size. Terrestrial 
gravity is eliminated from Equation (1) by immersing 
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the “planet” and its ferrofluid “ocean” in an immisci- 
ble silicone oil/Freon mixture of the same density. 
Thus the earth gravity is removed from the dynamics 
of the ferrofluid/oil interface and the only dynamically 
active force there is the radial magnetic gravity. The 
entire apparatus can rotate, and waves are forced on the 
ferrofluid surface by exterior magnets. 



Figure 2. Side view schematic cross section 
of the cylindrically symmetric experiment. A 
cylindrical tank holds a silicone oil and Freon 
mixture (2) surrounding a water-based ferrofluid 
layer (1) that is held by the magnetic body force 
onto the surface of a PVC ball covering a cylin- 
drical permanent magnet. p\ = P 2 to eliminate 
terrestrial gravity. The entire apparatus rotates 
at angular frequency, Q. 

Figure 3 is a perspective view photograph of the 
apparatus. The external disk-shaped magnet at left 
forces a wave with maxima symmetrically displaced off 
the equator and a local minimum on the equator. 
Waves with a single maximum on the equator can be 
forced by inverting the forcing magnet polarity. A 
Plexiglas wall just visible on the right side in the 
black ferrofluid serves as a continental boundary for the 
ferrofluid ocean for studies of wave reflection discussed 
below. The entire ferrofluid ball is about 35cm in di- 
ameter, twice the size of reference [3], and the surround- 
ing oil/Freon mixture has a viscosity 1/5 as large, giv- 
ing more than an order of magnitude reduction in vis- 
cous effects. Careful temperature control of the water in 
the outer square box is required to maintain uniform 
density between the ferrofluid and oil/freon mixture as 
the expansion coefficients of the two fluids are different 


The biggest improvement in technique over [3] is 
in the wave visualization. Fluorescing dye is added to 
the oil/Freon mixture and an argon ion laser generates 
the horizontal light sheet visible in Figure 3. This 
sheet can be scanned vertically. Viewed from above, 
the experiment is a black circle with wave deformations 
surrounded by a light background. A contour of the 
image intensity at any light sheet position gives the 
surface of the ferrofluid “ocean” at that “latitude”. Ra- 
dial displacements of the waves as a function of longi- 
tude are obtained by subtracting the contour line posi- 
tions from a no-motion contour at that laser sheet lati- 
tude. The experiments are run by traversing the forc- 
ing magnet with the laser sheet height fixed and im- 
ages are frame grabbed to obtain a time series at one 
latitude. The experiment is then re-run with another 
laser-sheet height to build up a full picture of the three- 
dimensional wave structure in the upper hemisphere of 
the ball as a function of time 



Figure 3. Photograph of the laboratory ex- 


periment. See text above for a description. 

RESULTS 

We concentrate here on results of laboratory studies 
of waves that are important in Earth's atmosphere and 
especially the ocean. In large scale atmospheric and 
oceanic fluid dynamics, only the local vertical projec- 
tion of the planetary rotation, / = 2f2sin(latitude), is 
important because the atmosphere and ocean are very 
thin compared to their horizontal extent and stratifica- 
tion effects further inhibit vertical motions. An impor- 
tant class of waves in geophysical fluid dynamics de- 
pends on the latitudinal variation of / [4, 5]. These 
Rossby waves are evident for example in the structure 
and dynamics of the atmospheric jet stream and Atlan- 
tic Gulf Stream. At the equator / changes sign, and a 
special class of these waves is trapped to equatorial 
latitudes but carry energy and momentum east and 
west [5]. These equatorial Kelvin and Rossby waves 
are an integral part of the dynamics of El Nifio. This 
ferrofluid experimental technique is the only laboratory 
realization of these equatorial waves. To get oceanic 
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scaling in the laboratory, the experiment must rotate 
rapidly (4-second rotation period) so that the wave 
speed is slow compared to the planetary rotation speed 
as in the ocean. 

In [3] observations of equatorial Kelvin and Rossby 
waves were compared to theoretical predictions [4, 5] 
and satellite observations [6, 7]. These experiments 
were conducted on the full sphere with visualization cf 
wave height only on the equator. Using the improved 
experiment with lower friction and better visualization, 
we have extended the results in [3] to verify the latitu- 
dinal shape of the westward propagating equatorial 
Rossby waves. These waves have their largest signal 
off the equator, something like the wave shown in Fig- 
ure 3. 

In the Pacific Ocean, eastward propagating Kelvin 
waves eventually run into the South American coast. 
Theory predicts that some of the wave energy should 
scatter into coastal -trapped Kelvin waves that propagate 
north and south along the coast [8]. Some of this 
coastal wave energy might then scatter into mid- 
latitude Rossby waves that propagate back westward. 
Satellite observations of the Pacific Ocean sea-surface 
temperature and height seem to show signatures of 
westward propagating Rossby waves at 20N to SON 
latitude, 5 to 10 years after the 1982-83 El Nirlo [9]. 
The observational data is difficult to interpret unambi- 
guously owing to the large range of motions that fill 
the ocean at shorter timescales. 



-200 -100 0 100 200 
Longitude West (deg) 


Figure 4. Wave height at 0°N latitude for 
non-rotating experiment. 



Figure 5. Wave height at 0°N latitude for 4- 
second rotation period experiment. 
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Figure 6. Wave height at 60°N latitude for 
4-second rotation period experiment. 
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This series of reflections giving eastward, north- 
ward, and then westward traveling waves is observed 
cleanly in the laboratory experiments, confirming the 
theoretical expectations. Figures 4, 5, and 6 are all 
contours of wave height vs. longitude and time fcr 
experiments with a vertical wall (continent) placed at 
longitude -130°W. In the experiments, a wave with 
maximum amplitude on the equator is forced by an 
external magnet moving in the co-rotating direction 
(eastward). The magnet is closest to the experiment at 
0°W longitude. Gaps in the laser-light sheet projec- 
tion system give rise to the contour noise from -20 °W 
to 10°W and 160°W to -170°W longitudes. In Fig- 
ure 4 the experiment is not rotating and the wave 
height on the equator is plotted. The wave reflects 
from the wall and propagates at the same speed back to 
the west. Much of the wave energy propagates toward 
the pole in the non-rotating case (not shown). In con- 
trast for the rapidly rotating case (Figures 5 and 6), the 
eastward wave propagates to the wall along the equator 
in Figure 5 and then disappears at about 40 seconds. 
Plots at higher latitude show the wave propagating 
northward until in Figure 6 at 60°N, the wave reaches 
maximum amplitude along the wall at 45 seconds and 
then propagates westward, more slowly than the east- 
ward wave. This westward propagation at speeds 
slower than the equatorial Kelvin wave speed fits the 
classic long Rossby wave signature. 

FUTURE WORK 

The linear and nonlinear characteristics of traveling 
or standing waves on a liquid surface are a classic prob- 
lem in fluid dynamics. Lateral boundaries have a 
strong effect on the spatial structures (fronts, pulses, 
solitons, etc.) observed in nonlinear propagating wave 
systems [10, for an extensive review]. It is of great 
theoretical interest to contrast wave behavior with and 
without lateral boundaries. With our spherical mag- 
netic body force, a straight circumferential channel 
formed by boundaries at ± latitudes is an experimental 
system for free-surface capillary-gravity waves that has 
periodic boundary conditions in one direction. Simi- 
larly, the entire spherical surface is boundary-less and 
periodic in two-dimensions. In an equatorial channel 
geometry, the contact lines are of the same length (as 
opposed to in an annular geometry) and this symmetry 
makes comparison to theory easier. We plan next to 
study various capillary-gravity waves using this appa- 
ratus which is the laboratory version of the standard 
theoretical periodic-boundary construct 

Standing capillary-gravity waves are excited if the 
gravitational force oscillates in time. This is normally 
accomplished by vibrating a fluid container vertically. 
These “Faraday waves” have become, along with 
Rayleigh-Benard convection, a canonical system for 
studies of non-equilibrium pattern formation, spatio- 
temporal evolution, and chaos [10]. Our spherical 
geometry, with parametric oscillation obtained by os- 


cillating the magnetic field, is a very clean system for 
investigating such parametrically excited waves, with- 
out complications of boundary shapes and contact an- 
gles which can determine patterns and mode competi- 
tion dynamics even very for from the boundaries [11]. 
We’ve designed an electromagnet version of the appa- 
ratus of Figures 2 and 3 to provide an oscillating 
“gravity”. We should have sufficient electromagnetic 
driving strength in the frequency range 0.01 to 10 Hz, 
to generate Faraday wave patterns by parametric excita- 
tion with up to a few tens of waves in a circumference 
on the sphere or somewhat fewer fitting around inside a 
polar basin. The spherical system has the periodic 
boundaries of the standard theoretical development [11] 
but which have not been previously possible experi- 
mentally. 

Continuing the geophysical applications, ferrofluid 
experiments offer the potential for conducting labora- 
tory studies of thermally driven oceanic flows over a 
substantial part of a sphere. These slow deep overturn- 
ing oceanic motions are highly constrained by the dy- 
namical influence of planetary vorticity advection. 
Such effects cannot be studied in the terrestrial labora- 
tory, when the all-important case of continuous thermal 
stratification is considered, because of the strong grav- 
ity induced cells that have no oceanic analogs. Previ- 
ous GFFC-type [12] experiments are also ill suited to 
oceanic flow modeling because they are too viscous 
(thin gaps). We will address the key technical issues 
in building a relatively simple but effective micrograv- 
ity based experiment to enable laboratory study of con- 
tinuously stratified spherical flows in complex basins 
spanning a latitude range of at least plus or minus 60 
degrees. We realize that such an experiment must be 
conducted in microgravity, and our terrestrial demon- 
stration experiment will be contaminated by terrestrial 
gravity induced flows. The main technical issues to 
explore are: 1) how to get a large enough magnetic 
buoyancy frequency in an apparatus of sufficient size, 2) 
how to stratify and force motion in the experiments, 
and 3) how to visualize the flows effectively. 
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